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We  use  a  theoretical  mesoscopic  model  accounting  for  memory  and  nonlocal  effects  in  thermoelectricity 
in  order  to  investigate  how  the  figure-of-merit  in  cylindrical  thermoelectric  nanodevices  is  conditioned 
in  frequency-dependent  situations.  Two  different  situations,  regarding  the  relative  values  of  the 
particles'  mean-free  path  and  the  characteristic  size  of  the  system,  are  analyzed.  It  is  shown  that  in 
both  situations  the  performances  of  the  thermoelectric  devices  are  reduced.  However,  nonlocal  effects 
may  be  used  as  an  aiding  tool  to  have  less  marked  reductions  in  those  performances. 
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1.  Introduction 

Alternative  power  sources  based  on  energy  harvesting  are 
promising  candidates  to  substitute  batteries  due  to  their  ability 
to  extract  unlimited  power  from  the  environment  or  secondary 
processes,  as  well  as  to  attain  fully  autonomous  systems  without 
periodical  human  interventions.  Due  to  the  large  amount  of 
residual  heat  yielding  from  the  current  energy  generation  tech¬ 
nology  based  on  fossil  fuels,  thermoelectric  energy  harvesters  have 
received  special  attention  in  recent  years.  Thermoelectric  devices 
offer  a  very  attractive  source  of  energy  since  they  do  not  have 
moving  parts,  create  pollution,  or  make  noise. 

In  practical  applications,  the  definition  of  a  “good  thermo¬ 
electric  device”  is  usually  related  to  the  dimensionless  product  ZT, 
with  T  being  the  operating  temperature,  and  Z  the  so-called  figure- 
of-merit,  defined  as 


wherein  e  is  the  Seebeck  coefficient,  cse  the  electrical  conductivity, 
and  X  the  total  thermal  conductivity. 

In  thermoelectric  materials  heat  is  carried  both  by  phonons  and 
by  electrons,  and  by  definition,  in  Eq.  (1)  the  total  thermal 
conductivity  is  such  that 

X  =  Xp  +  Xe  (2) 

with  Xp  being  the  phonon  contribution  to  the  thermal  conductivity, 
whereas  Xe  means  the  electron  contribution  to  it  [1-3]. 
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Since  the  higher  ZT,  the  higher  the  efficiency  of  a  thermo¬ 
electric  device,  in  order  to  widen  the  applications  of  thermo¬ 
electric  power  generators,  in  the  last  decades  there  has  been  a 
tremendous  amount  of  researches  to  improve  the  values  of  ZT 
beyond  those  of  bulk  materials,  which  show,  instead,  low  efficien¬ 
cies.  However,  ZT  has  remained  approximately  equal  to  1  for  the 
past  several  decades  in  the  case  of  archetype  materials  at  all 
temperature  ranges.  These  materials  include  antimony  (Sb2Te3) 
and  bismuth  tellurides  (Be2Te3)  for  room  temperature  applica¬ 
tions,  lead  telluride  (PbTe)  at  moderate  temperatures,  and  silicon- 
germanium  (SiGe)  alloys  at  high  temperatures. 

One  of  the  primary  challenges  in  developing  advanced  thermo¬ 
electric  materials  is  decoupling  e,  oe  and  X  which  are  typically 
strongly  interdependent  in  such  a  way  that  an  increase  in  e  usually 
results  in  a  decrease  in  oe,  and  a  decrease  in  oe  produces  a  decrease 
in  the  electronic  contribution  to  X,  following  from  the  Wiede- 
mann-Franz  law. 

Indeed,  if  the  characteristic  dimension  of  the  material  (or  the 
system)  is  shortened,  the  new  variable  of  length  scale  also 
becomes  available  for  the  control  of  the  materials'  properties.  In 
particular,  as  the  system  size  decreases  and  approaches  nanometer 
length  scales,  new  opportunities  are  allowed  to  vary  the  afore¬ 
mentioned  parameters  quasi-independently.  Nanomaterials, 
therefore,  provide  an  interesting  avenue  to  achieve  this  goal,  for 
example,  making  nanocomposites,  adding  nanoparticles  to  a  bulk 
material,  or  employing  one-dimensional  nanostructures  [4,5]. 
Nanosystems  offer  the  possibility  of  an  additional  control  of  the 
transport  coefficients  [6-9],  For  instance,  whenever  the  character¬ 
istic  size  of  the  system  is  comparable  to  the  mean-free  path  (mfp) 
f  of  the  different  heat  carriers  (phonons,  electrons,  holes,  etc.)  it  is 
known  that  a  thermal  conductivity  reduction  can  be  realized  over 


A  Sellitto  /  Physica  B  456  (2015)  57-65 


a  wide  temperature  range,  or  the  power  factor  can  be  increased  at 
the  same  time  by  increasing  e  more  than  ae  is  decreased  [10,11]. 

Although  from  experimental  evidences  it  is  very  clear  the 
importance  of  using  nanotechnologies  in  thermoelectricity,  the 
search  of  a  very  good  thermoelectric  device  is  still  far  from  its  final 
solution.  This  principally  because  the  physics  at  nanoscale  still 
presents  several  dark  points,  as  for  instance  the  role  played  by 
memory,  nonlocal  and  nonlinear  effects  [12,13,6,8,7], 

In  particular  memory  effects  may  drastically  influence,  for 
example,  the  behavior  of  nanodevices  at  high  frequencies 
[14,15],  In  fact,  in  the  simplest  description  of  relaxational  effects 
in  the  bulk  for  phonon  heat  transport,  the  heat  flux  q  is  given  by 
the  Maxwell-Cattaneo  equation 

rpDtq+q=  — 2VT  (3) 


where  r  is  the  relaxation  time  due  to  the  resistive  interactions 
(momentum  not  conserved)  between  the  heat  carriers,  and  the 
symbol  Dt  means  the  material  derivative.  The  Fourier  transform  of 
this  equation  leads  to  an  effective  frequency-dependent  thermal 
conductivity  of  the  form 


(4) 


which  points  out  that  the  higher  the  frequency  co,  the  smaller  the 
thermal  conductivity. 

This  result  may  be  interesting  for  thermoelectric  applications, 
since  small  values  of  the  thermal  conductivity  could  lead  to  an 
enhancement  of  Z,  according  to  Eq.  (1). 

Indeed,  the  analysis  and  modelization  of  heat  transport  in 
nanosystems  are  more  complex  since  it  is  also  required  to  pay  a 
special  attention  to  the  boundary  conditions  [3],  The  interest  on  wall 
effects,  related  to  phonon- wall  interactions  [16,17],  recently  surged 
with  the  synthesis  of  nanowires  with  rough  walls  [18,19],  Further¬ 
more,  in  thermoelectricity,  one  has  to  monitor  also  the  frequency 
behavior  of  other  material  functions.  Therefore,  in  the  present  paper, 
from  a  phenomenological  point  of  view,  we  explore  the  frequency 
dependence  of  Z  in  cylindrical  nanowires  in  order  to  point  out 
whether  the  performances  of  thermoelectric  devices  may  be 
enhanced  by  coupling  relaxational  and  nonlocal  effects  in  high- 
frequencies  situations,  or  not.  It  is  important  to  note  that  the 
approach  used  in  the  present  paper  rests  on  a  mesoscopic  level,  as 
our  main  aim  is  to  provide  a  first  rough  (but  simple  to  obtain) 
estimation  of  the  response  of  a  thermoelectric  device  in  frequency- 
dependent  situations.  In  this  sense,  the  present  analysis  should  be 
only  viewed  as  a  first  step  towards  more  detailed  microscopic 
scrutinies. 

The  structure  of  the  paper  is  the  following.  In  Section  2  we 
introduce  transport  equations  with  relaxational  and  nonlocal 
effects  in  heat  and  electric  transport.  In  Section  3  we  derive  the 
frequency-dependent  behavior  of  the  figure-of-merit  by  modeling 
the  interactions  of  the  different  heat  carriers  with  the  boundaries. 
In  Section  4  we  draw  the  main  conclusions. 


2.  Enhanced  thermoelectric  equations  with  relaxational  and 
nonlocal  effects 

On  microscopic  grounds  both  electrons  and  phonons  may  be 
viewed  as  a  free-particle  gas  in  a  box  [20],  In  nonequilibrium 
mechanics,  the  statistical  behavior  of  a  thermodynamic  system  far 
from  its  thermodynamic  equilibrium  is  described  through  the 
Boltzmann  transport  equation  (BTE),  which  in  the  relaxation¬ 
time  approximation  reads: 

dtf+  ■  ■  Vv/=  (5) 


wherein  the  subscripts  r  and  v  in  the  nabla  operators  represent 
the  variables  of  the  gradient,  i.e.,  they  are  the  position  and  the 
velocity  of  a  set  of  particles,  respectively.  Moreover,  in  Eq.  (5) 
/(r;v;t)  is  the  probably-density  function,  f0  represents  the  equili¬ 
brium  distribution  of  the  carriers,  and  frr;  k)  is  their  relaxation 
time,  k  being  the  wave  vector  of  the  particle.  In  the  BTE,/0  is  given 
by  the  Bose-Einstein  distribution  function  (BEdf)  in  the  case  of 
phonons,  i.e., 

/o  =  (eWV_i)-1 

wherein  h  =  h/( 2ri)  with  h  being  the  Planck  constant,  v  is  the 
angular  frequency,  and  kB  is  the  Boltzmann  constant.  In  the  case  of 
electrons,  instead,  f0  is  expressed  by  the  Fermi-Dirac  distribution 
function  (FDdf),  i.e., 


with  ei  being  the  energy  of  the  single-particle  state,  and  fie  is  the 
chemical  potential. 

Indeed,  it  is  possible  to  find  several  situations,  as  for  example 
whenever  the  quantities  (e, ■—//<,)  and  ftp  are  much  larger  than  kBT, 
in  which  one  can  ignore  the  +  1  in  the  denominator  of/0,  in  order 
that  the  BEdf  and  the  FDdf  reduce  to  the  Boltzmann  distribution 
function  [20],  In  these  cases,  on  intuitively  ground,  the  solution  of 
the  BTE  both  for  phonons,  and  for  electrons  would  lead  to 
equations  showing  the  same  mathematical  behavior. 

Starting  from  these  considerations  and  in  accordance  with  the 
basic  principles  of  Extended  Irreversible  Thermodynamics  [7,6],  in 
Refs.  [21,3]  the  following  generalized  transport  equations  to 
describe  heat  and  electric  transport  with  thermoelectric  coupling 
have  been  considered: 

TPDtq(p>  +  q®  =  -2pVT+«?2(V2q®+2VV  •  q®)  (6a) 

reDtq(e)  +  q(e)  =  -(4+ e!7<7e)VT +  <?2(V2q®  +  2VV  •  q(e)) 

+IIaeE  (6b) 

reDti+i  =  cre(E  — eVT)  +  if2(V2i+2VV  ■  i)  (6c) 

In  these  equations  i  is  the  electric-current  density,  E  is  the 
electric  field,  and  n  is  the  Peltier  coefficient.  The  basic  idea  lying 
under  these  equations  is  that  the  local  heat  flux  shows  both  a 
phonon  contribution  q®,  and  an  electron  contribution  q(c|  in  such 
a  way  that  q  =  q®  +  q(e). 

Referring  the  readers  to  Section  4  for  more  comments  about 
Eqs.  (6),  here  let  us  only  comment  about  the  different  material 
functions  included  therein.  In  Eqs.  (6)  rp  and  te  represent,  respec¬ 
tively,  the  relaxation  time  due  to  phonons  interactions  and  that 
due  to  electrons  interactions.  In  more  detail,  the  relaxation  time  zp 
may  be  related  to  the  resistive  mechanisms  between  the  different 
particles  in  such  a  way  that,  through  the  Matthiessen  rule,  it  is 
given  as 

*p  1  =  AT  1  + Tf  1  +  Td- 1  +  T~\  +  T-Jp 

where  ru  is  the  relaxation  time  of  umklapp-phonon  collisions, 
Tj  the  relaxation  time  of  phonon-impurity  collisions,  id  the  relaxa¬ 
tion  time  of  phonon-defect  collisions,  and  re_p  is  the  electron- 
phonon  scattering-time  [22],  Similarly,  the  electron  relaxation 
time  re  may  be  defined  as 


with  Te_e  being  the  electron-electron  scattering-time  [23,24],  In 
principle,  in  defining  both  tp  and  re  further  relaxation  times,  account¬ 
ing  for  the  phonon-wall  and  electron-wall  interactions,  should  be 
taken  into  account.  In  our  approach  these  interactions  will  be 
introduced  by  means  of  suitable  boundary  conditions  in  Section  3. 

Moreover,  in  Eq.  (6a)  €p  means  the  mfp  of  phonons,  which  is 
related  both  to  resistive,  and  to  the  normal  scattering  of  phonons 
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[25]  since  C^  =  vjt„ rp/5,  t„  being  the  relaxation  time  of  normal 
collisions,  and  vp  being  the  average  phonon  speed.  In  Eqs.  (6b)  and 
(6c),  instead,  €e  is  the  mfp  of  electrons  in  the  bulk  [26],  which  is 
dominated  by  electron  and  phonon  scattering  in  such  a  way  that 
fe  =  vF re,  with  vF  as  the  Fermi  velocity.  Indeed,  €e  can  be  also 
expressed  by  the  Drude  formula  [27], 

In  Eqs.  (6)  the  second-order  spatial  derivatives  of  the  fluxes 
allow  to  account  for  nonlocal  effects,  whereas  their  time  deriva¬ 
tives  introduce  memory  and  relaxational  effects.  It  is  easy  to 
recognize  that  in  the  absence  of  nonlocal  and  memory  effects, 
from  Eqs.  (6)  the  classical  thermoelectric  equations 
q=  -2VT+77i 
i  =  <7e(E  —  eVT) 

are  recovered  once  Eqs.  (6a)  and  (6b)  are  summed. 

For  the  sake  of  a  theoretical  completeness,  since  the  model 
Eqs.  (6)  lie  on  the  idea  that  the  heat  carriers  form  a  gas-like 
mixture  flowing  through  the  crystal  lattice,  we  observe  that  each 
heat  carrier  obeys  the  same  balance  equations  as  the  single 
component  of  the  mixture  [3],  namely,  the  following  further 
equations  hold: 

Dtu<P>  +  V  -q<P>  =  0  (8a) 

Dtu(e>  +  V  ■  qfe|  =  E  ■  i  (8b) 

where  u(m  and  u(e)  are  the  partial  internal  energies  (of  phonons 
and  electrons,  respectively),  and  the  term  E  •  i  accounts  for  the 
Joule  dissipation.  The  rate  of  the  electric  charge  density  c><e>, 
instead,  obeys  the  charge  conservation  law,  i.e., 

Dt<?(e)  +  vi  =  0  (9) 

which,  together  with  Eqs.  (8),  determines  the  evolution  of  u®,  u(e) 
and  Q(e>,  once  the  evolutions  of  the  corresponding  fluxes  q(p),  q(e) 
and  i  are  known. 

In  the  next  section,  we  will  analyze  how  the  dynamical  terms  in 
Eqs.  (6)  (i.e.,  TpDtq(p>,  TeDtq(e>,  and  reDti),  complemented  with 
suitable  boundary  conditions  (accounting  for  the  interactions 
between  the  different  heat  carriers  and  the  walls),  may  influence 
the  figure-of-merit  in  frequency-dependent  situations. 


3.  Frequency  dependent  figure-of-merit  in  cylindrical 
nanowires 

Owing  to  the  fast  progress  in  the  fabrication  of  nanoscale 
devices  and  the  predictions  of  an  improvement  in  the  figure-of- 
merit  of  low  dimensional  semiconducting  structures,  many  efforts 
have  been  made  in  the  thermoelectric  field  via  engineering  at  the 
nanoscale.  Although  the  main  research  line  rests  on  the  search  of 
suitable  new  materials  with  high  values  of  Z,  also  new  strategies 
should  be  desirable.  Motivated  by  this,  in  Ref.  [3]  the  authors  have 
studied  the  size  dependency  of  the  figure-of-merit,  complement¬ 
ing  Eqs.  (6)  with  suitable  boundary  conditions  accounting  for  the 
interactions  between  the  different  heat  carriers  and  the  lateral 
walls  of  the  system  in  steady-state  situations.  That  analysis  clearly 
pointed  out  the  role  played  by  nonlocal  effects,  and  their  possible 
use  to  enhance  the  performances  of  a  thermoelectric  device. 
However,  at  nanoscale,  also  memory  and  relaxational  effects  may 
be  important,  especially  in  unsteady-state  situations  [14],  There¬ 
fore,  in  the  present  section,  we  try  to  estimate  how  the  figure-of- 
merit  in  cylindrical  nanowires  may  depend  on  the  frequency  of 
perturbations.  In  particular,  in  order  to  analyze  the  possible 
consequences  of  accounting  for  relaxational  effects  in  practical 
applications,  here  we  apply  the  theoretical  model  in  Eqs.  (6)  in  the 
concrete  situation  of  a  nanowire  which  is  characterized  both  by  a 


(constant)  circular  transversal  section  of  radius  R,  and  by  a  long¬ 
itudinal  length  L  much  larger  than  the  latter.  We  will  perform  our 
analysis  in  the  following  two  possible  cases: 

(i)  when  R  is  such  that  fe<R<fp\ 

(ii)  when  R  is  such  a  way  that  R  <  t?e,  €p. 

In  both  cases  we  assume  that  the  temperature  difference 
(applied  through  the  ends  of  the  system)  and  the  electric  field 
may  vary  in  the  following  sinusoidal  way: 

AT  (or,  t)  =  AT0e’"'t  (10a) 

E(ar,  t)  =  E0ei'"c  (10b) 

w  being  the  angular  frequency  of  the  perturbations,  while  AT0  and 
E0  are  stationary  reference  levels.  Eqs.  (10)  imply  that  both  the 
temperature,  and  the  electric  field  are  homogeneous  across  every 
transversal  section  of  the  nanowire,  but  with  an  amplitude 
changing  in  time.  Moreover,  the  former  depends  on  the  long¬ 
itudinal  position,  whereas  the  latter  is  independent  of  it. 

3.1.  Case  of  study:  i?e<R<Cp 

In  Si  or  Bi2Te3  the  phonon  mfp  is  generally  one  order  of 
magnitude  larger  than  the  electron  one  [28,19].  Therefore,  when 
these  materials  are  used  in  building  up  thermoelectric  devices,  a 
realistic  assumption  is  to  assume  that  the  radius  of  the  transversal 
section  is  larger  than  fe,  but  smaller  than  fv. 

In  this  case  the  phonons  undergo  the  hydrodynamic  regime, 
and  the  electrons  the  usual  resistive  regime.  In  hydrodynamic 
regime  the  transport  is  ballistic  [14],  and  in  the  evolution  equation 
of  the  corresponding  heat  flux,  i.e.,  in  Eqs.  (6a),  the  second-order 
spatial  derivatives  play  the  main  role,  in  such  a  way  that  their 
contribution  is  predominant  with  respect  to  that  of  the  heat  flux 
itself  [17,3],  Moreover,  in  such  a  situation,  the  heat-flux  profile  in  a 
given  transversal  section  is  strongly  dependent  on  the  character¬ 
istic  size  of  the  system  [29,17],  In  the  resistive  regime,  instead,  the 
transport  is  diffusive  [14],  and  the  second-order  spatial  derivatives 
of  the  diffusive  flux  can  be  neglected,  since  the  main  role  in  the 
evolution  equation  is  played  by  the  diffusive  flux  itself.  In  this  case, 
the  dependence  of  the  single  diffusive  flux  on  the  characteristic 
size  is  realistically  negligible. 

From  the  practical  point  of  view,  previous  considerations  allow 
to  conclude  that  whenever  fe<R<fp  Eqs.  (6)  may  be  reduced  as 


xpDtqfp)  =  -ipvr+^p(V2q(p)+2VV  ■  q(p)) 

(11a) 

reDtqm  +  q,e)  =  -(4  +  eIIoe)VT+noeE 

(lib) 

teDti+i  =  <re(E— eVT) 

(11c) 

In  order  to  keep  the  calculations  down  to  a  reasonable  level,  in 
using  Eqs.  (11)  we  assume  that  both  the  heat  flux,  and  the 
electrical  current  may  only  flow  longitudinally  through  the  sys¬ 
tem.  Moreover,  in  the  present  case,  we  assume  that  the  bulk 
profiles  for  the  different  diffusive  fluxes  are  given  by 

q$\r,co;t)  =  qw(*1^  je* 

(12a) 

qf(®;t)  =  q^e'w 

(12b) 

tb  (®;  t)  =  ib,0e''”t 

(12c) 

that  is,  we  assume  that  they  have  amplitudes  which  periodically 
deviate  in  time  from  the  stationary  reference  levels  q®,  q®  and 
zb(| ,  but  they  are  independent  of  the  longitudinal  position.  It  is 
worth  noticing  that  Eq.  (12a)  prescribes  that,  in  a  given  section  and 
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time,  q{f}  essentially  keeps  the  parabolic  behavior  corresponding  to 
the  Poiseuille  phonon  flow  [29,17],  whereas  q®  and  ib  are  homo¬ 
geneous.  These  assumptions  follow  from  the  considerations  above 
about  the  different  transport  regimes. 

We  also  note  that  qjp) ,  q£°  and  ibo  in  Eqs.  (12)  are  related  to  A T0 
and  £0  in  Eqs.  (10).  In  fact,  combining  Eqs.  (10)  and  (12),  by 
straightforward  calculations,  we  have 


+  2ja)C?Vp 

1+jWp(R4^2) 


AT0 


(13a) 


!'b,o  = 


(e  ad  (eE)  ne+enae\AT0  (  TIoe  \_ 

-<  +<  -  (T+jRr)TT+(T+JS;)E« 

-  ’b,0  +  'b.0  ~~  (i  +j£BTe)  i-  +  (l  +jarte)  E° 


(13b) 

(13c) 


wherein  c®  is  the  specific  heat  at  constant  volume  of  phonons, 
which  may  be  related  to  the  temperature  through  the  usual 
constitutive  assumption  u®  =  c®T.  It  seems  worth  observing  that 
in  the  whole  generality,  u®  should  also  be  also  q  and  i  dependent. 
However,  although  more  refined  assumptions  may  be  introduced, 
here  we  use  that  assumption  in  order  to  deal  only  with  the 
simplest  situation,  emphasizing  the  physical  idea.  In  particular, 
that  constitutive  relation,  combined  with  Eq.  (8a),  has  been  used 
to  express  VV  •  q®  in  terms  of  the  temperature  gradient  in 
deriving  Eq.  (13a). 

Since  the  electrons  may  be  dragged  both  by  a  temperature 
gradient,  and  by  an  electric  field,  in  Eq.  (13b)  the  two  different 
contributions  to  the  bulk  electron  heat  flux  qj®  have  been 
explicitly  expressed:  that  only  due  to  the  temperature  gradient 
(i.e.,  q^AT> )  and  that  only  due  to  the  electric  field  (i.e.,  q^). 
Analogously,  in  Eq.  (13c)  appear  the  two  different  contributions  to 
the  bulk  electric-current  density  ibi0,  namely,  i^AT)  which  is  due  to 
the  temperature  gradient,  and  t®  which  is  due  to  the  presence  of 
the  electric  field. 

As  we  are  analyzing  a  situation  in  which  the  phonons  undergo 
the  hydrodynamic  regime,  we  have  to  pay  a  special  attention  on 
the  interactions  between  them  and  the  lateral  walls  [29,30],  This 
can  be  done,  for  example,  by  assuming  that  the  full-local  phonon 
contribution  to  the  heat  flux  is  given  as 


q®  =  q®  +  q® 


(14) 


phonon  Knudsen  layer,  in  fact,  the  electron-wall  interactions 
should  be  felt  in  a  further  thin  layer  near  the  walls,  the  character¬ 
istic  size  of  which  in  of  the  order  of  €e.  This  zone,  which  for 
homogeneity  of  nomenclature  we  call  electron  Knudsen  layer,  is 
less  wide  with  respect  to  the  phonon  Knudsen  layer,  since  we  are 
assuming  that  ee  <  R.  Therefore,  in  this  case,  we  do  not  commit  a 
significant  error  in  admitting  that  the  wall  contributions  to  q®  and 
i  are  vanishingly  small. 

These  considerations  allow  also  to  understand  that  the  particle- 
wall  interactions  are  strictly  related  to  nonlocal  effects,  since  they 
do  not  have  any  relevance  whenever  nonlocal  effects  may  be 
neglected. 

In  view  of  the  dynamical  form  of  q®,  in  the  same  way  as  above 
we  assume  that  the  wall  contribution  is  given  by 

q®(a>;t)  =  q®  e**  (15) 

wherein,  in  a  first-order  approximation  in  space,  the  reference 
amplitude  value  q®0  can  be  obtained  by  the  following  constitutive 
equation  [29,17,30]: 

r ®Dtq®  +  q®  =  Cptp  (16) 

In  Eq.  (16)  t®  means  the  relaxation  time  accounting  for  the 
phonon-wall  collisions.  In  principle,  it  may  be  estimated,  for 
example,  by  the  usual  Matthiessen  rule  as 

'w^'spJc  +  'd® 


where  Tspec  and  rdiff  are,  respectively,  the  characteristic  times  of 
specular  and  diffusive  reflections  [29],  Indeed,  as  a  function  of  the 
characteristic  features  of  a  wall  (i.e.,  the  square  value  of  the 
roughness  fluctuations  and  the  average  distance  between  rough¬ 
ness  peaks  [31]),  in  Ref.  [32]  the  authors  provided  a  possible  way 
of  estimating  that  relaxation  time  in  terms  of  the  total  frequency 
of  collisions  between  phonons  and  walls. 

Moreover,  in  Eq.  (16)  Cp  means  a  dimensionless  coefficient 
accounting  for  the  rate  of  reflections  (specular  and  diffusive)  of 
phonons  [29,3],  In  Ref.  [29]  the  authors  provided  a  simple  model 
to  estimate  it  as  a  function  of  the  temperature  and  of  the  wall 
features.  In  Ref.  [33],  instead,  the  following  more  refined  model  for 
Cp  may  be  found: 


Cp 


2  [(3  -jf3)  3  (1  -/2)1 
'3[  I  2~ Knp~J 


where  q®  is  the  wall  contribution  which  is  progressively  tuned 
down  far  from  the  walls.  Note  that  in  general,  in  accounting  for 
boundary  effects,  it  should  be  admitted  that  q®  takes  the  value  q[pl 
in  the  bulk,  and  the  value  q^  at  the  walls.  However,  since 
boundary  conditions  strongly  determine  the  behavior  of  a  system, 
one  could  nevertheless  expect  that  the  thermal  conductivity  is 
only  dependent  on  q®  in  the  bulk,  and  on  the  value  of  q®  at  the 
boundary  of  the  wall.  Indeed,  from  the  physical  point  of  view,  the 
boundary  is  never  a  true  surface,  but  rather  a  thin  layer  (the  so- 
called  Knudsen  layer),  whose  dimension  is  of  the  order  of 
magnitude  of  the  mfp  of  the  heat  carriers.  On  the  other  hand, 
since  R  <  «?p,  the  Knudsen  layer  pervades  the  whole  system,  so  that 
the  heat  flux  on  the  boundary  (i.e.,  in  the  phonon  Knudsen  layer)  is 
superposed  to  the  heat  flux  in  the  bulk.  This  fact  suggests  to 
assume  Eq.  (14)  to  estimate  the  overall  heat  flux,  and  suppose  that 
the  thermal  conductivity  depends  on  it.  Moreover,  whenever 
R  <  €p,  we  are  also  allowed  to  assume  that  q®  attains  the  same 
value  in  each  point  of  the  transversal  section. 

No  wall  contributions  are  introduced,  instead,  for  the  full-local 
values  of  q®  and  i  which,  therefore,  coincide  with  their  corre¬ 
sponding  bulk  values.  Along  with  previous  observations  about  the 


with  /  =  min[l  /Knp,  1],  Knp  =  t?p/R  being  the  phonon  Knudsen 
number,  and  £  a  suitable  accommodation  parameter.  This  coefficient 
represents  the  portion  of  the  total  wall-colliding  particles  that  are 
diffusively  reflected  back  by  wall  and  have  a  bulk  velocity  equal  to 
the  wall  velocity  after  collisions.  In  principle,  Cp  values  should  be 
taken  in  the  range  [0;1  ].  In  fact,  since  Cp  also  turns  out  the  relative 
importance  of  q®  with  respect  to  q®  in  Eq.  (14),  in  general  the  Cp 
values  should  be  small  in  such  a  way  that  q{f’  always  remains  the 
main  contribution  to  the  full-local  phonon  heat  flux  in  Eq.  (14). 

Note  that  Eq.  (16)  has  to  be  meant  as  a  relaxational  general¬ 
ization  of  the  boundary  condition  used  in  Ref.  [3]  (see  Eqs.  (11) 
therein)  in  the  spirit  of  the  Maxwell-Cattaneo  equation  (3), 
wherein  a  relaxational  term  (including  the  time  derivative  of  the 
heat  flux)  is  added  to  the  usual  Fourier  heat-transport  equation. 
Though  in  the  usual  applications  the  term  r®Dtq®  is  vanishingly 
small,  for  high  frequencies  it  may  be  obviously  important,  and 
therefore  its  incorporation  in  Eq.  (16)  seems  logical.  Referring  the 
readers  to  Section  4  for  further  comments  about  Eq.  (16),  we 
observe  that  a  simple  and  systematic  derivation  of  that  equation 
can  be  found  in  Ref.  [30]  by  using  the  tools  of  Extended  Irrever¬ 
sible  Thermodynamics  [7], 
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The  coupling  of  Eqs.  (13a)  and  (15)  with  Eq.  (16)  turns  out  that 
the  amplitude  of  the  wall  contribution  to  the  phonon  heat  flux  is 

c,  %+3-tVjVr. 


'v’°  21<np  y  i  +ja)Tw  J  L 
If  we  define  the  following  effective  quantities 


(4  |  cllo,.ff  -  Rei 


L  /Qg.AT) 


I  AT 


1} 


(18b) 

(18c) 


then  all  the  considerations  above  allow  to  obtain  the  following 
effective  coefficients: 

1"-2J'(S)2{ln[1+fe)  1  -'"““'"““(si) 


(19a) 

4  =l  +  (®re)2 

(19b) 

jefrfJL  'Vj 
e  1  +  (arref 

(19c) 

e^  =  e 

(19d) 

ncff  =  n 

(19e) 

wherein  ap  =  r(f/rp.  For  the  sake  of  formal  simplicity,  in  Eq.  (19a) 
we  expressed  the  phonon  relaxation  time  as  tp  =  > fp/vp,  and  we 
used  the  usual  Ziman  limit  for  the  phonon  bulk  thermal  con¬ 
ductivity,  i.e.,  Xp  =  c®TpV^/ 3. 

It  seems  interesting  to  note  that  the  theoretical  model 
described  by  Eqs.  (11)  prescribes  that  only  the  thermal  and  the 
electrical  conductivities  show  effective  values  that  differ  from 
their  corresponding  bulk  values.  Moreover,  as  expected  whenever 
ife  <R  <  Xp,  among  them  only  the  phonon  contribution  to  the 
effective  thermal  conductivity  shows  the  influence  of  the  radius. 

It  is  easy  to  see  that  in  the  low-frequency  limit  (that  is, 
whenever  anp  -*  0  and  a)te->0)  Eqs.  (19a)-(19c)  reduce  to  the 
effective  values  obtained  in  Ref.  [3]  whenever  the  radius  of  the 
transversal  section  is  in  the  range  [fe;  ifp\,  that  is, 


In  Fig.  1  we  plot  the  frequency-dependent  behaviors  of  xf1,  xf* 
and  nff  (referred  to  their  corresponding  bulk  values)  arising  from 
Eqs.  (19a)-(19c)  in  the  case  of  a  p-doped  sample  of  Bi2Te3  at  300  K. 


Furthermore,  the  combination  of  Eqs.  (1)  and  (19a)-(19c),  is  used 
to  plot  in  the  same  figure  also  the  behavior  of  the  effective 
frequency-dependent 

7eff_^W 

referred  to  its  corresponding  bulk  value. 

The  values  of  the  different  material  functions,  related  to 
phonons  and  electrons,  we  used  in  the  computation  are  summar¬ 
ized  in  Table  1.  For  the  sake  of  computation,  we  assumed  R=  2  nm. 
From  those  values  we  estimated  the  following  relaxation  times: 

tp=^  =  3.5x10_13s 

re  =  —  =  2.8  x  10-33  s 
vF 


This  means  that  whenever  onp  ->  1 ,  then  co  is  of  the  order  of 
THz.  Therefore,  in  Fig.  1  co  spans  the  usual  range  of  admissible 
frequencies.  Since  we  are  not  aware  of  experimental  data  to 
estimate  the  parameter  accounting  for  the  phonon-wall  reflec¬ 
tions,  we  used  two  different  admissible  values  for  it,  that  is, 
Cp=0.1  and  Cp=0.5. 

As  it  can  be  seen  from  Fig.  1,  the  higher  o>,  the  smaller  xf  and 
of1.  These  results  are  a  natural  consequence  of  the  electron 
transport-regime  type  (i.e.,  the  resistive  regime).  In  the  case  of 
the  phonon  transport  regime,  instead,  the  relevance  of  nonlocal 
effects  leads  to  an  enhancement  of  xff  for  increasing  co.  The 
combination  of  these  behaviors  leads  to  an  effective  figure-of- 
merit  which  is  substantially  decreasing  for  increasing  frequency, 
whatever  the  value  of  Cp  is.  However,  the  smaller  Cp,  the  less 
marked  the  decreasing  in  Z.  Therefore,  the  theoretical  model  in 
Eqs.  (6)  suggests  that  the  performances  of  a  thermoelectric  device 
tend  to  go  worse  in  frequency-dependent  situations. 

3.2.  Case  of  study:  R  <  te,  €v 

Due  to  the  incessant  developments  in  nanotechnology,  it 
should  be  also  possible  to  face  with  situations  in  which  the 
characteristic  size  of  the  device  is  smaller  both  than  the  phonon 
mfp,  and  than  the  electron  one.  In  these  situations  all  the  heat 
carriers  undergo  the  hydrodynamic  regime,  and  Eqs.  (6)  become 

TpDtq®  =  -  /lpVT+^p(V2q,p)  +  2VV  ■  q(p))  (20a) 

TeDtq<e)  =  -(Xe  +  cnoe)VT+  *f2(V2q(e)  +  2VV  •  qw)  +  /7oeE  (20b) 
TeDti  =  cre(E  —  eVT)  +  if2(V2i+2VV  •  i)  (20c) 

In  this  case,  in  order  to  study  the  frequency-dependent  effec¬ 
tive  figure-of-merit  we  can  still  assume  that  qb  is  given  by 
Eq.  (12a),  which  substantially  leads  to  the  frequency-dependent 
effective  phonon  thermal  conductivity  (19a),  once  the  phonon- 
wall  interactions  are  introduced  by  a  wall  contribution  as  in 
Eq.  (16). 

Different  treatments,  instead,  are  needed  for  q\f  and  ib.  In  fact, 
since  in  this  case  also  the  electron  transport  is  ballistic,  we  can  no 
longer  assume  that  those  fluxes  show  profiles  which  are  homo¬ 
geneous  in  each  transversal  section  for  a  fixed  time.  Therefore, 
along  with  previous  considerations  about  Eq.  (12a),  here  we 
assume  that  qf*  and  ib  have  the  following  profiles,  respectively: 


(21a) 
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Fig.l.  Frequency-dependent  behaviors  of  if1,  >epK,  nf  and  of  Z*,fr  in  the  case  ofap-doped  sample  of  Bi2Te3  at  T=300  K:  theoretical  results  arising  from  the  theoretical-model 
equations  (11). 


Table  1 

Values  of  the  material  functions  for  a  p- doped  sample  of  Bi2Te3  at  300 1<  [28,19], 


4  4  vr  Vp  4  (m)  4  (m) 

(Wm-'r1)  (Wm-’K-1)  (ms-1)  (ms-1) 


2.4  1.6  3.26  xlO3  8.43  x  103  9.1  xlO10  3.0  x  10“: 9 


ib(®;  r;  t)  =  ib,0  ^  (21b) 

namely,  we  are  assuming  that  and  ib  also  keep  the  parabolic 
behavior  corresponding  to  the  Poiseuille  flow  [29,17], 

In  Eqs.  (21 )  the  reference-level  amplitudes  and  ib  0  are  still 
strictly  related  to  the  temperature  difference  and  the  electric  field 
through  the  end  of  the  system,  but  in  different  forms  with  respect 
to  those  predicted  by  Eqs.  (13b),  and  (13c).  In  fact,  if  we  combine 
Eqs.  (20b)  and  (20c)  with  Eqs.  (21),  by  direct  calculations  we  have 
that 


(4  +  elide)  +  2/<uCve)  (f g 

AT0 

Tide 

,  .  (r2-v*\ 

L  + 

1+jWe(^) 

(22a) 


g  e 

AT0 

G 

L  1 

(22b) 


wherein  c®  is  the  specific  heat  at  constant  volume  of  electrons, 
which  we  relate  to  the  average  temperature  T  through  the 
constitutive  assumption  u(e)  =  Cy}T  for  a  formal  simplicity. 
For  the  sake  of  clarity  in  reading  the  present  paper,  we  observe 
that  in  obtaining  Eqs.  (22)  we  combined  that  constitutive  assump¬ 
tion  with  Eqs.  (8b)  and  (9)  in  order  to  write 

VV  •  q(e)  =  V(  — C^DtT+E  •  i )^df)DtAT/L 

once  the  further  simplifying  assumption  of  vanishing  value  of 
Dtq(e)  in  Eq.  (9)  is  made. 

Since  now  we  are  assuming  that  R  is  also  smaller  than  t?e, 
here  we  can  no  longer  neglect  the  influence  of  the  electron- 
wall  interactions  [3],  In  other  words,  in  the  present  case  we 
have  that  also  the  electron  Knudsen  layer  pervades  the  whole 
transversal  section.  Therefore,  in  analogy  with  Eqs.  (14), 
here  we  assume  that  the  full-local  values  of  q®  and  ib  are 
given  by 

qW  =  q<f  +  qW  (23a) 

i  -  ib  +  iw  (23b) 

To  cope  with  high-frequencies  situations,  still  in  analogy  with 
what  supposed  for  the  phonon-wall  heat  flux,  we  assume  that  the 
wall  contributions  q®  and  iw  in  Eqs.  (23)  show  the  following 
dynamical  behaviors: 

q®(®;  t)  =  q„\  e!°*  (24a) 

iw(ca;  t)  = 


(24b) 
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with  the  reference  amplitudes  q^]0  and  iw>0  satisfying  the  dynami¬ 
cal  constitutive  assumptions 


T($Dtq$+q$  =  Ce/fe 


(25a) 


r®Dtiw 


(25b) 


wherein  is  a  relaxation-time  due  to  electron-wall  interactions,  and 
Ce  turns  out  the  rate  of  reflected  electrons.  They  substantially  play 
analogous  roles  of  and  Cp,  respectively.  In  particular,  once  Eqs.  (24) 
are  introduced  into  Eqs.  (25),  by  straightforward  calculations  one  has 


,w-“  2oKne[^ 


+  enae  +  2jmcifVA  A  T0 
1  +jan$  )  L 


V  +J®Tw)/ 


(26a) 


Ce  1 
+  4Kne|l  ~  (aeoney 


1} 


whereas  the  effective  Seebeck  and  Peltier  coefficients  still  keep 
their  corresponding  bulk  values.  Note  that  in  deriving  Eqs.  (27)  we 
expressed  the  electron  relaxation  time  as  te  =  €e/Vf,  and  the 
electron  bulk  thermal  conductivity  as  Ae  =  tf)Tev}/3.  Moreover,  in 
those  equations  we  set  ae  =  r<$  /ve. 

It  is  easy  to  see  that  in  the  low-frequency  limit  (that  is, 
whenever  o>rp -* 0  and  a>Te-> 0)  Eqs.  (27)  reduce  to  the  effective 
values  obtained  in  Ref.  [3]  whenever  R  <  ee,  ifp,  that  is, 


-2Kne[(i+j ,W$)  L°+ 


(26b) 


where  Kne  =  ife/R  is  the  electron  Knudsen  number. 

Then,  by  means  of  Eqs.  (18b)-(18e),  whenever  R  <  €e  we  obtain 
that  the  frequency-dependent  effective  values  of  the  electron 
thermal  conductivity  and  of  the  electrical  conductivity  in  Eqs. 
(19b)  and  (19c)  become 


-  12fflre  arctan 


51} 


(27a) 


In  Fig.  2  we  plot  the  frequency-dependent  behaviors  of  ,  /pf 
and  (referred  to  their  corresponding  bulk  values)  arising  from 
the  theoretical  model  in  Eqs.  (6)  in  the  case  of  a  p-doped  sample  of 
Bi2Te3  at  300  K,  whenever  the  radius  of  the  transversal  section  is 
smaller  both  than  €v,  and  In  that  figure  the  behavior  of  the 
effective  frequency-dependent  Zefr,  referred  to  its  corresponding 
bulk  value  Z,  is  also  plotted.  For  the  sake  of  computation,  in  this 
case  we  assumed  R=0.5  nm  which  leads  to  the  following  relaxa¬ 
tion  times  of  heat  carriers-walls  interactions: 

tJS?  =  ^  =  5.9  x  10-14  s 

— sstl.5  X  10~13  s 
vF 


Fig.  2.  Frequency-dependent  behaviors  of  I'f ,  Ajf ,  and  of  Zeff  in  the  case  of  a  p-doped  sample  ofBi2Te3  at  T=300K:  theoretical  results  arising  from  the  theoretical-model 

equations  (20). 
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whereas  zp  and  re  get  the  same  values  as  in  Section  3.1.  Moreover, 
since  we  are  not  aware  of  experimental  data,  in  this  case  we  still 
used  two  different  values  for  the  parameter  accounting  for  the 
phonon-wall  reflections  (i.e.,  Cp=0.1  and  Cp=0.5),  and  two  differ¬ 
ent  values  for  the  parameter  accounting  for  the  electron-wall 
reflections  (i.e.,  Ce=0.1  and  Ce=0.5). 

As  it  can  be  seen  from  Fig.  2,  increasing  values  of  a>  lead  to 
increasing  values  in  both  the  thermal  conductivity  2pff  and  /.ff , 
whereas  the  effective  electrical  conductivity  of1  shows  a  slight 
decreasing  behavior.  We  explicitly  note  that  the  opposite  beha¬ 
viors  of  and  aff  one  finds  in  Fig.  2  are  not  against  the 
Wiedemann-Franz  law  since  at  nanoscale,  as  we  observed  in 
Section  1,  the  characteristic-length  variable  allows  additional 
control  on  the  different  material  functions.  The  combination  of 
these  behaviors  leads  to  an  effective  figure-of-merit  which  is 
substantially  decreasing  for  increasing  frequency,  whatever  the 
values  of  Cp  and  Ce  are.  However,  the  smaller  the  parameters 
accounting  for  the  reflection  of  phonons  and  electrons  at  the  walls, 
the  less  marked  the  decreasing  in  Z.  Therefore,  the  theoretical 
model  in  Eqs.  (6)  suggests  that  the  performances  of  a  thermo¬ 
electric  device  tend  to  get  worse  in  frequency-dependent  situa¬ 
tions  even  when  the  radius  of  the  transversal  section  is  smaller 
both  than  the  phonon  mfp,  and  than  the  electron  one. 


4.  Conclusions 

Thermoelectrics  is  a  meaningful  energy-conversion  technology, 
which  could  generate  electrical  power  from  a  temperature  differ¬ 
ence  and  also  could  produce  active  cooling  by  inputting  direct 
electrical  current.  Indeed,  the  low  energy-converting  efficiency  by 
thermoelectric  system  has  been  the  bottleneck,  which  limited  the 
applications  to  niche  markets  where  the  reliability  and  simplicity 
are  more  important  than  the  performance.  The  revitalization  of 
the  thermoelectric  industry  has  been  due  to  nanotechnologies,  the 
use  of  which  requires  a  deeper  scientific  understanding  of  the 
physics  at  nanoscale.  From  the  theoretical  point  of  view,  it  is  well- 
known,  for  example,  that  resistive,  capacitive  and/or  inductive 
properties  of  nanomaterials  show  interesting  time-dependent 
features  when  subject  to  time-dependent  perturbations,  but  how 
they  may  influence  the  performances  of  thermoelectric  devices  are 
still  an  open  question. 

In  the  present  paper,  starting  from  a  thermodynamic  model  of 
nonlocal  transport  equations  for  the  phonon  and  the  electron 
contributions  to  the  heat  flux,  and  for  the  electron  contribution  to 
the  electric  current,  including  thermoelectric  effects,  by  means  of 
the  figure-of-merit  we  have  analyzed  how  time-dependent  per¬ 
turbations  may  influence  the  performances  of  a  thermoelectric 
device. 

Since  Eqs.  (6)  represent  the  basic  tool  of  our  theory,  it  seems 
advisable  to  spend  more  words  about  them. 

Eq.  (6a)  is  very  well-known  in  nonequilibrium  thermody¬ 
namics  as  the  Guyer-Krumhansl  (G-K)  equation  [25],  and  follows 
from  the  solution  of  the  linearized  BTE  for  phonons  [34,7,35],  It 
reduces  to  the  classical  Fourier  law  whenever  the  terms  rpDtq^ 
and  ifp(V2q<p)  +  2VV  •  q®)  are  neglected  therein.  Along  with  the 
observations  at  the  very  beginning  of  Section  2,  in  the  present 
paper  we  assumed  that  the  electron  contribution  to  the  local 
heat  flux  is  ruled  by  Eq.  (6b),  with  additional  terms  accounting 
for  the  Peltier  effect.  Whenever  the  terms  zeDt q(e)  and 
if2(V2 q(e)  +  2VV  •  q(e))  can  be  ignored  in  Eq.  (6b),  one  recovers  the 
classical  equation  describing  the  Peltier  effect  [36,6],  namely, 
q(e)  =  -  (2e  +  cJlae)  V T  +  Tlae E.  Finally,  since  the  electric  current  is 
produced  by  the  electric  charges  moving  through  the  crystal 
lattice,  we  also  assumed  that  their  evolution  is  governed 
by  a  G-K  equation  type,  that  is,  Eq.  (6c).  That  equation  reduces 


to  i  =  ce(E— eVT)  whenever  the  relaxational  term  r,Dti  and  the 
second-order  nonlocal  contribution  ^(V2i+2VV  •  i)  are 
negligible  in  it. 

The  results  contained  in  Ref.  [37]  about  the  compatibility  of 
Eqs.  (6)  with  second  law  of  thermodynamics  allow  to  claim  that 
those  equations  do  not  violate  any  admissible  physical  principle. 

In  order  to  account  for  the  interactions  between  heat  and 
electric  carriers  and  the  lateral  walls,  in  our  theoretical  model  we 
complemented  Eqs.  (6)  with  suitable  boundary  conditions  describ¬ 
ing  the  wall  contributions  to  the  different  diffusive  fluxes.  The 
establishment  of  boundary  conditions  in  the  presence  of  wall  slip 
has  been  a  subject  of  increasing  interest  in  fluid  mechanics, 
rheology  and  heat  transport  processes.  In  microfluidics,  for  exam¬ 
ple,  along  the  lateral  surfaces  of  a  cylindrical  duct  it  is  generally 
supposed  a  slip-flow  condition  (i.e.,  the  fluid  velocity  does  not 
vanish  along  it),  whereas  the  behavior  of  the  fluid  in  the  center  of 
the  duct  is  governed  by  the  usual  Navier-Stokes  equations.  In  a 
first-order  approximation,  the  velocity  vw  of  the  fluid  on  the  wall 
of  the  microchannel  is  generally  estimated  as 


wherein  vb  represents  the  velocity  of  the  fluid  in  the  bulk,  f  is  the 
mfp  of  the  fluid  particles,  c  is  a  positive  constant  related  to  the 
properties  of  the  walls,  and  f  means  the  outward  normal  direction 
to  the  boundary  [38-40],  Along  with  the  idea  that  the  heat/electric 
carriers  form  a  gas-like  collection,  replacing  the  fluid  velocity  with 
the  corresponding  dissipative  flux  J  of  the  heat/electric  carriers,  in 
the  case  of  heat  transport  at  nanoscale  from  Eq.  (28)  one  may  infer 
that  J  on  the  walls  is  given  by 

Jw  =  (29) 

ot, 

It  is  easy  to  recognize  that  Eqs.  (16)  and  (25)  reduce  to  Eq.  (29) 
in  steady-state  situations,  or  whenever  the  relaxational  effects  due 
to  the  interactions  between  the  heat/electric  carriers  and  the  walls 
are  vanishing.  A  rigorous  theoretical  study  about  the  compatibility 
of  Eqs.  (16)  and  (25)  may  be  performed  as  in  Ref.  [30]  by  imposing 
that  the  entropy  production  is  non-negative  not  only  in  the  bulk, 
but  also  in  the  bounding  surface. 

In  the  case  of  a  cylindrical  nanodevice,  we  have  analyzed  two 
different  situations  regarding  the  relation  between  the  radius  of 
the  transversal  section  and  the  mfps  of  the  different  heat  carriers. 
In  both  situations  our  theoretical  model  (6)  predicts  that  Z  is 
reduced  for  increasing  frequency  of  perturbations.  Indeed,  our 
analysis  suggests  that  nonlocal  effects  may  be  used  to  bring  the 
reduction  in  Z  down  in  a  frequency-dependent  situation.  In  fact,  if 
in  the  synthesis  and  processing  of  a  nanomaterial  one  is  able  to 
address  the  different  particles  at  the  walls  in  such  a  way  that  only 
a  scant  number  of  them  suffer  reflections  (that  is,  if  Cp  and/or  Ce 
get  small  values),  than  the  reduction  in  the  performance  of  a 
thermoelectric  device  should  be  attenuated  in  frequency- 
dependent  situations.  From  this  perspective,  the  phenomenon  of 
backscattering  [41  ]  (i.e.,  the  backward  reflections  of  hitting  parti¬ 
cles  at  the  walls),  which  is  a  huge  problem  for  heat  conduction, 
may  be  used  as  an  aiding  tool. 
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